Modulational instability in periodic quadratic nonlinear materials 
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We investigate the modulational instability of plane waves in quadratic nonlinear materials 
with linear and nonlinear quasi-phase-matching gratings. Exact Floquet calculations, confirmed 
by numerical simulations, show that the periodicity can drastically alter the gain spectrum but 
never completely removes the instability. The low-frequency part of the gain spectrum is accurately 
predicted by an averaged theory and disappears for certain gratings. The high-frequency part is 
related to the inherent gain of the homogeneous non-phase-matched material and is a consistent 
spectral feature. 



The modulational instability (MI) of waveforms is 
a fundamental phenomenon in nonlinear media and 
is closely associated with the concept of self-localized 
waves, or solitons. Like solitons, MI occurs due to an 
interplay between nonlinear and dispersive or diffractive 
effects. It involves the exponential growth of weak per- 
turbations through the amplification of sideband frequen- 
cies, which causes a destructive modulation of the car- 
rier wave. Many diverse physical systems exhibit MI, 
e.g. fluids jl], plasmas nonhnear optics ||], molec- 
ular chains \%, and Fermi-resonant interface waves and 
waveguide arrays Q . In nonlinear optics MI may appear 
as a transverse instability that breaks up a broad optical 
beam, thereby acting as a precursor for the formation 
of stable bright spatial solitons [||. Conversely, the sta- 
ble propagation of dark solitons relies on the stability of 
the constant-intensity background and thus requires the 
absence of MI [§ . 

Here we study optical media with a purely quadratic 
(x*^^)) nonlinearity. These materials are of significant 
technological interest due their strong and fast cascaded 
nonlinearities which can support stable bright soli- 
tons in all dimensions Moreover, they provide the 
most general means of studying nonlinear processes, be- 
cause varying the phase mismatch between the funda- 
mental and second-harmonic (SH) waves changes the 
nonlinearity from being distinctly quadratic to effectively 
cubic. The generality is further enhanced when the prop- 
erties of the medium are modulated with long-period 
quasi-phase-matching (QPM) gratings. This periodicity 
not only allows one to tune the effective mismatch, but it 
also induces effective cubic nonlinearities [ p^pl| ], which 
may be engineered to different strengths and signs [Q. 

Except in special cases, MI is unavoidable in x'^^ mate- 
rials. Optical pulses can be modulationally stable only if 
the dispersion has opposite sign at the fundamental and 
SH frequencies p^ . Transverse inhomogeneities, as in 
X^^) waveguide arrays, may also eliminate MI How- 
ever, CW beams propagating in pure homogeneous x'^"* 



materials are always unstable 0|. Naturally, a cubic 
(Kerr) nonlinearity is always present in x'^^ materials, 
and if defocusing and sufficiently strong, it may elimi- 
nate MI Unfortunately, the cubic nonlinearity in 
conventional x*-^-* materials is usually focusing, though it 
may be a factor in the recent observation of apparently 
stable quadratic optical vortex solitons [|l5| . 

We consider transverse MI of CW beams in x'^' mate- 
rials whose linear and nonlinear properties are both mod- 
ulated and make a detailed study of the instability gain 
spectra. Because such spectra can now be measured in 
the laboratory studying their profiles has direct ex- 
perimental relevance. Our purpose is to identify regimes 
with no MI, or with a maximum gain that is small enough 
to allow spatial dark solitons to propagate stably over 
experimentally relevant distances. It is known that, in a 
system with cubic nonlinearity, periodic dispersion IT^ ] 
(as in dispersion-managed fibers) or nonlinearity can 
reduce the growth rate of fluctuations. We find that, in 
quadratic materials, the periodicities also drastically al- 
ter the MI gain spectrum but they do not entirely remove 
the instability. 

We numerically find exact plane- wave solutions and ap- 
ply Floquet theory to determine their MI gain spectrum. 
Numerical simulations, with exact plane-wave solutions 
seeded with noise as initial conditions, confirm the Flo- 
quet gain spectra. In the regime of efhcient quasi-phase- 
matching, we find that the MI gain spectrum contains 
two distinct and well-separated features with fundamen- 
tally different physical origins. The low-frequency part of 
the gain, which to our knowledge has not been predicted 
before, is accurately predicted by an averaged theory and 
disappears for certain grating modulations. The high- 
frequency part of the spectrum is related to the inher- 
ent gain of the non- phase- matched material (i.e., with 
no gratings) and appears to be unavoidable. However, 
because they are consistently small, the high-frequency 
peaks can be ignored under a less stringent definition of 
experimental stability. This gives the possibility of stable 
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plane waves as predicted by the averaged theory. 

The system under consideration consists of a CW beam 
(carrier frequency w) and its SH propagating in a loss- 
less ID slab waveguide under conditions for QPM type-I 
second-harmonic generation (SHG). Both the linear and 
quadratic nonlinear susceptibilities are periodic along the 
Z-direction of propagation, as illustrated in Fig. but 
are constant along the transverse X-direction. The mod- 
ulation of the refractive index is assumed to be weak: 
Anj(Z)/nj^l, where nj{Z) = hj + Anj{Z) and j refers 
to the frequency jto. For many materials and etching 
techniques, such a linear grating appears inevitably when 
the nonlinear susceptibility is altered We consider 
gratings for forward QPM only, such that the grating pe- 
riod is much longer than the optical period and we can 
neglect Bragg reflections. The evolution of the slowly 
varying beam envelopes is then given by |po| 

idzEi + \d\E^ + OLr{Z)E^ + x{Z)ElE2e^^ = 0, 

idzE2 + \d\E^ + 1a^{Z)E2 + x{Z)EIq-''^^ = 0, (1) 

where Ei=Ei{X,Z) and E2=E2{X, Z) are the envelope 
functions of the fundamental and SH, respectively. The 
coordinates X and Z are in units of the input beam 
width Xq and the diffraction length Ld—kiX^, respec- 
tively. The parameter P—AkL^ is proportional to the 
mismatch Ak=k2 — 2ki, where kj=jujnj/c is the average 
wave number. Thus f3 is positive for normal dispersion 
and negative for anomalous dispersion. The normalized 
refractive- index grating is aj{Z)=LdUjAnj{Z)/c, and the 
normalized nonlinear grating is xiZ)=EdLL!dcg{Z)/ (nic), 
where dcff =x''^V2 is given in MKS units. 

To consider a plane wave with longitudinal wave num- 
ber offset A and transverse wave number fl, we transform 
to the new variables w{x, z)^Ei {X, Z) exp {i^X — iAZ), 
v{x,z) = E2{X, Z) exp{2inx - 2iAZ + if3Z), x=^{X + 
nZ), and z=r]Z, where 7] = \A + for A ^ -^Q^ 

and 77=1 otherwise. We have also introduced a residual 
effective mismatch /3=/3 — kg. This gives the equations 

idzW + ^dlw -rw + a'^{z)w + x' {z)w* ve^'^'' = 

id^v + ^dlv ~av + 2oi^{z)v + x {z)w^ e-^''"'' = 0, (2) 



where r — sgn(A -I- ^i^^), cr = 
Note that r^O when A^-if^^ 



2r — P/i], and k — kg/rj. 
. Wc Fourier-expand the 



rescaled gratings, a^(z)=aj(Z)/77 and x'C-^) = x{Z)/r]: 
= 5"e , X = do + rf 2^ 5„e , (3) 

n n 

where, for the square waveform that we consider (Fig. ^), 
gn—'2s/{inn) for n odd and .g„=0 for n even. The sign 
s=sgn(K) is included to ensure that the expansions give 
the same grating form for both signs of k. The gratings 



will drive a periodic variation in z of spatial frequency k 
in the fields of system (||). We explicitly include this fast 
variation by expanding the fields in Fourier series also: 



(4) 



For plane-wave solutions {'w,v)—{ws,Vs), all the coeffi- 
cients w„ and Vn are constants determined by 

- {nn + r)wn + a[ Qn-iwi + Dn+i-p^iW^Vp = 0, 
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i,p 



{uK + Cr)Vn + 2a'2 ^ gn-m + ^ Dn-l-p+lWlWp = 

I 



l,p 



where Dn=d'gn + d^Sn, Sn being Kronecker's delta. We 
find the solution numerically using a standard relaxation 
technique based on Newton's method. 

We determine the linear stability of the plane-wave so- 
lutions by using exact Floquet theory [||. For perturba- 
tions of transverse frequency offset v. 

w{x, z) = Ws{z) + (5wi(z)e"*'''' + 5wl{z)e'-'''' , 

v{x, z) = v,{z) + 5vi{z)e-"''' + 5vl{z)e"'^ , (5) 

the linear growth is governed by the linearized equation 
9zP = M{z)'P, where P(2;) and the matrix M(z) are 
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and where we have defined the components a — 
r + a[{z), b = Vs{z)x{z)e'"^'' , c = Ws{z)x{z)e 
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2a'2{z). 



Because of the periodicity in a^-, x'j and Vs, the 
matrix M is also periodic with period Zp = 2i: /\k,\. To 
determine the stability, wc need only find the eigenval- 
ues Xi of the mapping 8T of the solution over one pe- 
riod: P(z -I- Zp) = STP{z). We calculate 6T numer- 
ically by integrating the linearized equation with a 
fourth-order Runge-Kutta method. If all the Xi lie on 
the unit circle for every v, then the solutions Ws and 
Vs are stable; otherwise there is a gain whose profile is 
g—max{di[ln{Xi) / Zp]} . We note that the actual growth 
rate at short distances may differ from g, depending on 
the overlap of the initial perturbation with the unstable 
eigenvector However, for larger distances (z > l/g), 
as we consider here, the growth rate will approach the 
limiting maximum value g. 

As a first example, we calculate the Floquet gain 
spectrum for a conventional domain-reversal grating in 
LiNbOs {a'i=a2=d'Q=0) for r=— 1, k= 100 and exact 
phase matching ((T=— 2). As Fig. ||(a) shows, the spec- 
tral features fall into two well-separated regions: two 
low-frequency bands around v=2.Q and v=2>A, and a 
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narrow high-frequency band at u=l4.2. For comparison 
Fig. ||(b) plots the gain for a GaAs/AlAs superstructure 
with a nonUnear grating etched though quantum- well dis- 
ordering ((ig/(i'=4.6) and a hnear grating chosen to 
be {a'l — a'2)/K—— 1/9.2. High-frequency gain bands are 
again present, but the low- frequency bands are absent. 
Further results show that such narrow, high-i^ bands are 
a general feature of the gain spectrum for QPM gratings. 
They are related to the inherent MI in the non-phase- 
matched, gratingless x^^^ material: Each gain peak that 
appears in the gratingless spectrum usually also appears 
in the Floquet spectrum, typically with the same height 
and spectral location. An exception is the gain from the 
symmetric grating case (dQ=0), as in LiNbOa, where the 
original 'homogeneous' peak is suppressed. Often, as in 
Fig. H(b) , this peak is split into two or more closely spaced 
components by the grating. 

The low-frequency features in the gain spectrum are 
accurately predicted by approximate averaged-field equa- 
tions for the DC components = (wp, vq). From 
Eqs. (Q) these equations are found to be |Tl|Jl0[| 



id. 



id. 



-1* 



- rw + pw*v + 7(|wp — = 0, 

av + p*w'^ + 2j\'w\'^v = 0, (7) 



where p = 2is{d' + 2dp{a'i — al^j k)It^ and 7 = (d!^ -\- 
— 8/7r^))/K. The averaged equations are valid for 
small residual mismatch (J5' <C /3) and small grating pe- 
riod (k 3> 1). They show that the linear grating may 
decrease the strength of the effective quadratic nonlin- 
earity, possibly even to zero |11 19]. More importantly, 



the nonlinear grating induces cubic nonlinearities. These 
can be defocusing (7 > 0) and may thus stabilize 
plane waves We find that Eqs. (|) indeed have mod- 
ulationally stable plane-wave solutions in a large regime 
for r = —1, as shown in Fig. §(a), where 7 — ^/\p\'^ 
gives the relative strength of the cubic nonlinearity. The 
figure also shows that, since the mismatch a cannot be 
large (otherwise the averaged equations are invalid), the 
effective quadratic nonlinearity \p\ must be reduced (i.e., 
7 increased) to get a stable solution. 

The case shown in Fig. ^(a) corresponds to having 
7 = 44.5, a = —2, which lies outside the region of stabil- 
ity predicted by the averaged equations. The gain profile, 
shown in Fig. |^(c), matches exactly the low- frequency 
part of the Floquet spectrum. The case in Fig. ^(b) cor- 
responds to p = 0, in which case the averaged equa- 
tions reduce to two coupled Schrodingcr equations. The 
plane- wave solution WQ~y^—r/^, vq—0 is stable for 7>0 
(r—~l), in agreement with the absence of low-frequency 
bands in the exact Floquet spectrum. 

Figure ^b) shows the maximum gain calculated from 
Floquet theory for a given modulation depth. Where 
the averaged equations predict MI, the Floquet gain is 
large, in agreement with the averaged results. However 



where the averaged equations predict stability, Floquet 
theory finds a residual gain, due to the high-frequency 
gain bands. However, because the residual gain is of- 
ten small and narrow, it may be neglected under a rea- 
sonable definition of experimental stability, which allows 
for a small growth over say 10 diffractions lengths (typ- 
ical length scale in experiments). Neglecting the high- 
frequency gain bands causes the results of Floquet and 
average-field calculations to coincide and whole regions in 
parameter space of "experimentally stable" plane-waves 
to appear. These regions are even larger than shown in 
Fig. ^ and also appear for r—+l and 7-=0. We demon- 
strate this experimental stability by simulating the prop- 
agation of dark solitons with Eq. (P, shown in Fig. ||(a). 
The parameters correspond to a nonlinear grating in- 
duced in GaAs/AlAs through quantum-well disordering 
|p2| , as in Fig. g(b) but with no linear modulation. Sta- 
ble propagation of the dark soliton can be seen for over 
60 diffraction lengths. Now because of the small gain in 
highly non-phase-matched materials, such stable propa- 
gation can also be achieved without the nonlinear grating 
(Fig. 1(b)). But the SHG efficiency is then very low, as 
evidenced by the relative SH intensity being a factor of 
10 less than with the grating. 

The Floquet spectral results are confirmed by numer- 
ical simulations of Eqs. (|l|). We launch plane waves 
seeded with Gaussian white noise to excite all spa- 
tial frequencies. The gain spectrum is then calculated 
from the Fourier-transformed fundamental w as g{i^) = 
{In\w{z2,i')\—ln\w{zi,i^)\)/{z2 — zi) As shown in 

Fig. H, the simulations agree with the Floquet curves, in 
both the \ow-iy (Fig. ||(b)) and high-i^ (Fig- ^(c)) regions. 

In summary, we have used an exact Floquet technique 
to find the MI gain spectrum in x^^^ materials with gen- 
eral QPM gratings. Due to the periodicity, novel high- 
and low-frequency gain bands appear in the spectrum. A 
simple approximate averaged-field theory correctly pre- 
dicts the low- frequency gain bands. The high-frequency 
peaks that are consistently present in the spectrum are 
typically small and may be neglected under a realistic 
definition of experimental stability. The predictions of 
the averaged and the exact Floquet analyses then agree, 
and stable plane waves and dark solitons become possi- 
ble, even under conditions close to phase matching. 

MI is yet to be investigated in x^^^ materials with both 
linear and nonlinear QPM gratings. However, it is now 
possible to directly probe the gain spectrum in experi- 
ments [ p^ , and thus the interesting spectral profiles pre- 
dicted here may be verified in the laboratory. 

This work is supported by the Danish Technical Re- 
search Council under grant no. 5600-00-0355. 
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FIG. 1. Normalized linear and quadratic nonlinear grat- 
ings, aj{Z) and x(^)i with domain length Lo=n /[kg\. 
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FIG. 2. (a) Region of MI predicted by the averaged-field 
theory with r=— 1 and p^O. (b) Maximum gain pre- 
dicted by Floquet theory for 5<k<30, with do/d'=5/3 and 
(a'l — a2)/K=0.158. The region of zero gain predicted by av- 
eraged-field theory is shaded. 
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FIG. 3. Floquet gain spectra for (a) the LiNbOs and 
(b) the GaAs/AlAs gratings, for a— — 2. The diamonds 
give the averaged-field results in (c) and the equivalent 
non-phase-matched homogeneous results in (d). 
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FIG. 4. Propagation of dark solitons with (a) a modulated 

nonlinearity [d'^/d' = 4.6, k = 112) and (b) no modulation 
(d' = 0). In both cases /3 = 100 and Ld ~ 200^im. 
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FIG. 5. (a) Gain profile calculated frompropagative simu- 
lations, for the same parameters as in Fig. ra(a). Comparisons 
with the Floquet theory (solid curve) are given in (b) and (c) . 
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